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Covariant R^ gauge fixing is notoriously difficult for large lattice volumes, large ^ and small 
Nc- We thoroughly test different convergence techniques, which allows the gauge fixing of lattice 
configurations with a total volume of (3.25 fm)^, up to = 0.5. We are able to study the gluon 
propagator in the infrared region and its dependence on the gauge fixing parameter ^. As ex¬ 
pected, the longitudinal gluon dressing functions stay constant at their tree-level value ^. Similar 
to the Landau gauge, the transverse R^ gauge gluon propagators saturate at a non-vanishing value 
in the deep infrared for all values of ^ studied. We compare with very recent continuum stud¬ 
ies and perform a simple analysis of the found saturation with a dynamically generated effective 
gluon mass. 
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Renormalized Gluon Propagator - |J. = 3 GeV 

P = 6.0 - combined 32^, 48^, 64*, 80* data 



Figure 1: (left) Gluon propagator and (right) gluon mass computed in the Landau gauge [3]. 


1. Introduction 

Dynamical mass generation in QCD and hadrons accounts for ~ 98% of the visible mass in 
the universe. Most of the mass of the light quarks is due to chiral symmetry breaking. An effective 
gluon mass (as the photon mass in a superconductor) is a possible signal of confinement [1, 2]. If 
it exists, the gluon mass should appear in the gluon propagator. However Green’s functions (prop¬ 
agators and vertices) depend on the gauge fixing. In fhe case of fhe Landau gauge, fhe fransverse 
part of fhe gluon propagafor safurafes in fhe IR, as illusfrafed in Fig. 1. This can be inferprefed as 
an evidence for a dynamically generafed mass [3]. While fhere is yef no unique way fo define fhe 
gluon mass, in general fhe gluon mass m^(p^) is a monofonically decreasing function, power-law 
suppressed in fhe UV. The quesfion is, whaf happens in ofher gauges? So far , we don’f really know. 

Here we sfudy in Laffice QCD fhe well known renormalizable-i§ covarianf gauges (including 
Landau gauge = 0 and Feynman gauge = 1 ). Reliable laffice calculafions in renormalizable- 
^ (R^) covarianf gauges have nof been sysfemafically pursued previously. gauge fixing (GF) 
laffice implemenfafion has in facf proven fo be quife complicafed [4]. Recenfly some success was 
achieved [5], however one sfill encounfers significanf convergence problems in realistic laffices, i e 
for larger GF pai'amefer ^ and laffice size, and smaller number of colours Nc and laffice coupling 

P- 

We report on our success [6] fo GF in a realisfic laffice and presenf fhe SU(3) gluon propagafor 
in R^ gauges for a relatively large laffice volume (3.25 fm)^ and a GF paramefer up fo = 0.5. 


2 . 

2.1 


gauges framework and gauge fixing algorithm 

R^ gauges 

In fhe confinuum, R^ GF is achieved by adding fo fhe SU(Ac) Yang-Mills action fhe ferm, 


5r,P= (Tx 


m\2 




( 2 . 1 ) 
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U config. 
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(a) Histogram of the number of iterations nec¬ 
essary to converge the R^ GF, with FFT steep¬ 
est descent minimization of the functional 9. 
We illustrate 200 uncorrelated conhgurations for 
^ = 0.3 and L = 32. After 50000 iterations we 
stop FFT and cycle through other GF techniques, 
OVR, STR, RGT until convergence is reached. 



(b) (top) The 32"^ values of V A"’ evaluated for a 
conhguration gauge fixed at ^ = 0.5, grouped in 
5000 bins, compared with the Gaussian A™ with 
standard deviation ~ 0.316. (Bottom) Plot 
of = V -A^ — A'^; the two distributions coincide 
within v/0 precision. 


Figure 2: Our approach to achieve the GF within the desired precision of 6 < 10 


where ^ is the GF parameter, b'” are Nakanishi-Lautrup multipliers and A™ = 
condition. Going on-shell, = A"*, and the GF action takes a Gaussian form, 

gauges are obtained when the linear condition is chosen 


A'"[A] is the GF 
( 2 . 2 ) 


A 


m 




(2.3) 


When the gluon propagator is decomposed in transverse / longitudinal components. 


Ajivitj) — {gpv ) + ) ) (2-4) 


Slavnov-Taylor identities ensure that q^Ai^ = ^ to all orders. 

In the lattice, we use the Wilson action and the gauge links Ufi, are related to the gauge fields. 


A^{x + e^l2) 


Up{x)-Ul{x) 

2/go 


traceless 


(2.5) 


In R^ gauges, besides the usual integration over the link variables {x) , one has to integrate over 
the fields A = where each A™ is a Gaussian distribution with variance ^ ^ . 

The procedure for GF requires fo gauge rofafe all link variables [4, 5], 


Up{x) ^ g{x)U^{x)g^{x + e^) 


(2.6) 


where g are elements of the S\J(Nc) gauge group. 

In the Landau gauge case, which is the —)■ 0 limit of the R^ gauges studied here, the GF is 
implemented minimizing the functional, —ReTr s{^)Up (.^)g^ i^ + ^p), since the minimization 
ofr any functional of A^ (v) directly leads to the condition V -A'” = 0. 
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(a) Gluon longitudinal dressing function (b) Gluon transverse R^ dressing function q^Aj, 
q^Ai^, compared to the theoretical value (renormalized at /r = 4.317 [GeV]). 


Figure 3: Guon propagator dressing functions. 


2.2 gauge fixing algorithm 

In the general case of a non-vanishing a more complicated functional than the one of Landau 
GF needs to be minimized. In practice, the gauge transformation, g = is built as a 

product of a sequence of infinitesimal gauge transformations. For each infinitesimal transformation, 
5gj = 1 -|- we minimize the functional with respect to w'". We directly utilize the 

linear gauge condition of Eq. (2.3), to define our functional 6, depending on U^{x) via Eq. (2.5), 


e = 


N,L^ 


£Tr[A(jc)A^(x)] , 


A^{x + e^l2)-A^{x-e^l2) -A{x) . 
4 


(2.7) 


When 0—5-0, then A{x) -5- 0 and we reach the desired gauge condition V-A"* = A’” of Eq. (2.3). 

Thus, choosing w'" = aA'”, where a is a relaxation pai'ameter to be optimized, should reduce 
A with a steepest descent method. Our goal is to converge to a vanishing A(x) in all lattice points 
X. Based on our experience of studying gluon propagators with Eandau GE, we aim at a very small 
Q < 10^'^. However, this turns out, in practice, to be very difficult in a realistic lattice! 

We opt to extend our fully parallel and very fast GPU codes for Landau GF [7], for a thorough 
optimization of all possible techniques to minimize 0. We combine three different GE techniques 
optimized in the Eandau case [7], 

• the Fast Fourier Transform - accelerated steepest descent (FFT), 


• Over Relaxation (OVR) 


• and Stochastic Relaxation (STR). 


Each one of theses techniques, with an optimized convergence parameter a has a similar R^ GE 
success rate of only ~ 75% for ^ = 0.3 and ~ 40% for ^ = 0.5. This is illustrated in Eig. 2a 

After a large number of tests performed in our GPU servers, we finally find a solution. Cycling 
through EET, OVR and STR, for our hardest case of = 0.5, we increase the convergence success 
rate up to ~90% . Einally, for the remaining 10% cases, we perform a random gauge transformation 
(RGT), and restart the combined algorithm, till convergence is total [6, 8], see Eig. 2b. 
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(a) We show he transverse propagator At (b) Ratio of R^ transverse propagators to the 
(renormalized at ji = 4.317 [GeV] ). Landau gauge transverse propagator. 


Figure 4: Transverse gluon propagators for different gauges. The grey crosses are computed 
with the Landau gauge for a volume of 80^ and provide an estimate for the volume effects expected 
at = 0 [11]. 


3. Results and comparison with continuum studies 

3.1 Results for the gluon propagator 

We now compute the gluon propagator for j3 = 6.0 in a 32^ lattice. Linear covariant gauges 
on the lattice requires, not only an integration over the ensemble of configurations U^{x) [9], but 
an additional integration over the — I Gaussian distributed A™ fields. For fhe A infegrafion, we 
consider 50 differenl A’s for each configuration U^. Then, for each configuration of U^{x) and 
A(x) field, GF is applied. The lattice gluon propagafor, a fwo-poinf correlafion funclion, reads 

(A™(^)A«(^')) = 5'™A^v(9)^'5(?+r) , (3.1) 

where we Fourier Iransform from fhe positions x fo fhe momenfa q. The fransverse and longifudinal 
SU(3) propagator form factors are [10, 11], 

^ ~ 9p9v/9^)(A^(^)A™(— g)), 

^ L 9m9v/9^(A"‘(?Mv (-?)>• (3-2) 

Analytically, since fhe longifudinal propagafor Al remains equal to fhe free level one. we 
should have < 7 ^Al = The values of ^ chosen are ^ =0.1, 0.2, 0.3, 0.4, and 0.5. Indeed a fif of fhe 
dafa fo a consfanf, show in Fig. 3a, yields ^ = 0. 103(2), 0.203(2), 0.302(3), 0.402(3) and 0.502(3) 
respectively. 

The R^ gluon fransverse dressing funcfion is fully dynamical and non-perfurbafive. The 
simulafed volume is (3.25 fm)^, large enough to resolve fhe onsef of fhe non-perfurbafive effecfs. 
For comparison, fhe Landau gauge resulfs obfained for a symmefric laffice of L = 80 and j3 = 6.0 
(gray crosses) are also plotted in Fig. 3b. 
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(a) Previous theoretical prediction of the gluon 
propagator [12]. 



(b) Reconstructed gluon mass in the gauges, 
using the Nielsen identities of Ref. [12]. 


Figure 5: Comparing with continuum studies and reconstructing the gluon mass. 


We plot the gluon transverse propagators in Fig. 4a. As in the Landau gauge, we find they 
show an inflection point, implying that the associated spectral density is not positive definite; it is 
interpreted as a manifestation of confinement. They have a marked tendency to flatten towards the 
small momentum region, thus providing strong evidence that also in the / 0 case the behaviour of 
the lower modes of the lattice gluon field are tamed by the dynamical generation of a (momentum- 
dependent) gluon mass. The Landau gauge data on very large volumes suggest that simulations on 
larger physical volumes suppress the gluon propagator in the infrared region. This would only lead 
to at most to a decrease of about ~ 10% of At at small momenta for the transverse propagator. 

We also plot the ratio, in Fig. 4b, of the transverse propagator to the Landau gauge propagator 
Aj~^ as a function of the momentum for the two values ^ =0.1 and ^ = 0.5. The data confirms an 
IR hierarchy such that At (slightly) decreases for increasing values of the gauge fixing parameter. 
The maximum difference is of ~ 10% for ^ = 0.5. 


3.2 Continuum studies, gluon mass 


Our data for the gluon propagator is comparable to the one predicted in Ref. [12], shown in 
Fig. 5a, who include Nielsen identities and a gluon mass. Applying the same analysis of of Ref. 
[12] to our R^ lattice propagators, we find that for small the R^ dynamical mass behaves like 










(3.3) 


where is the Landau gauge gluon mass and to lowest order in the gauge fixing parameter, 
a{^) = l+ai(§,c((§) = Cnii§. 

Using as input our R^=0 (Landau gauge), L = 32, j8 = 6.0 data for the gluon and ghost prop¬ 
agators we obtain the Landau gauge dynamical mass We then solve the renormalization 

group improved equation and determine as well as the numerical value Cm ~ 0.32. We also deter¬ 
mine the numerical value of a i by requiring that the dynamical mass equals the value of A^r ^ (0) for 
the corresponding value of We obtain cq « 0.26. Finally including the statistical errors on the 
lattice propagators, we obtain the gluon running mass for the different R^, plotted in Fig. 5b. 
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4. Conclusions and outlook 

From the numerical point of view, our most intensive task [6, 8], is the gauge fixing due to 
the large number of GF’s required and algorithmic issues. Up to large ^’s and volumes, with a 
proper combination of various steepest descent methods, we solve the GF in gauges . We also 
compute the lattice SU(3) gluon propagators in gauges, for a lattice volume large enough to 
access the IR dynamics. Our Aj{q^) propagators are ~ similar for ^ = 0.0 — 0.5 : an inflection 
point in the few hundreds MeV region and a saturation in the IR. Our propagators are in agreement 
with very recent continuum analytic studies, and we estimate the dynamically generated R^ mass. 
This analysis suggests that dynamical gluon mass generation is a common feature of all R^ gauges 
in SU(3) Yang-Mills theories. 
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